Large deviations of ergodic counting processes: a statistical mechanics approach 
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The large-deviation method allows to characterize an ergodic counting process in terms of a 
thermodynamic frame where a free energy function determines the asymptotic non-stationary sta- 
tistical properties of its fluctuations. Here, we study this formalism through a statistical mechanics 
approach, i.e., with an auxiliary counting process that maximizes an entropy function associated 
to the thermodynamic potential. We show that the realizations of this auxiliary process can be 
obtained after applying a conditional measurement scheme to the original ones, providing is this 
way an alternative measurement interpretation of the thermodynamic approach. General results 
are obtained for renewal counting processes, i.e., those where the time intervals between consecutive 
events are independent and defined by a unique waiting time distribution. The underlying statistical 
mechanics is controlled by the same waiting time distribution, rescaled by an exponential decay mea- 
sured by the free energy function. A scale invariance, shift closure, and intermittence phenomena 
are obtained and interpreted in this context. Similar conclusions apply for non-renewal processes 
when the memory between successive events is induced by a stochastic waiting time distribution. 

PACS numbers: 05.70.Ln, 05.40.-a, 42.50.Lc, 02.50.-r 
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I. INTRODUCTION 

In different scientific disciplines the measurement real- 
izations of a given process consist of a set of random 
points distributed along the real number line. These 
"stochastic point processes" [ll-Q can be characterized 
through different statistical objects and techniques. One 
of the most usual is to count the (random) number of 
events up to a given point. For example, one can count 
up to a given time the number of photons emitted by a 
quantum optical system 0413 j the number of electrons 
transported through a nanoscopic structure (T^ - [l6| . or 
the number of spikes produced by a neuronal system 
[T7| - |2]| . just to name a few. From these examples it 
becomes clear that the underlying (non-equilibrium) dy- 
namics that lead to the random processes may be of very 
different nature and complexity. 

A counting process can statistically be characterized 
by a set of functions, each one defining the probability 
of n-events occurring up to a given time. These objects 
completely characterize the fluctuations (moments and 
cumulants) of the measurement realizations. The behav- 
ior (time-dependence) of the counting probabilities is not 
universal and depends on each specific case. Neverthe- 
less, after a characteristic transient time that depends 
on each system, a general characterization is expected to 
apply. In fact, when the probabilities of a system scales 
as exponential functions (asymptotic regime), the large 
deviation (LD) approach [22| allows to read the scaling 
rates in the context of an equilibrium thermodynamics 
formalism. In this frame [22- 29] , the entire (asymptotic) 
statistical properties of a process can be obtained from a 
free energy function j30 | and characterized, in particular, 
from its associated thermodynamic response functions 
(dynamical order parameters), i.e., its derivatives with 
respect to an intensive parameter (conjugate field). The 



thermodynamic frame describes in a unified way both 
the typical fluctuations of a system, such as those that 
can be fit in a Gaussian description [central limit theo- 
rem (CLT)], but also large deviation fluctuations that go 
beyond the previous regime. 

The LD theory and its associated thermodynamic 
frame were recently applied for characterizing photon 
counting processes in different quantum optical systems 
[28l [29} . The thermodynamic responses develop proper- 
ties such as scale invariant points, phase transitions and 
finite-size effects, which in turn indicate drastic changes 
in the statistical properties of the counting process. 

While previous results [2^, [2^ confirm that the LD 
method provides a deep description of a stochastic count- 
ing process, its associated thermodynamic structure can 
only be analyzed after specifying a given problem or sit- 
uation. Furthermore, while the full thermodynamic for- 
malism depends on the conjugate (counting) field it is 
not completely clear how this dependence may become 
observable or explicitly defined from the statistics of the 
counting process. The main goal of this paper is to shed 
light on this issue and to demonstrate that some general 
features of the thermodynamic frame can be character- 
ized without knowing the specific properties and under- 
lying dynamics (classical, quantum, phenomenological) 
of the counting process. The present analysis applies to 
ergodic (unidirectional) processes. 

Instead of focusing on the thermodynamic properties, 
which depend on each particular problem, here we ana- 
lyze its associated statistical mechanics, i.e., an auxiliary 
counting process that maximizes an entropy function re- 
lated to the free energy of the thermodynamic frame. We 
show that its realizations can be obtained after applying a 
conditional measurement scheme (Fig. 1) to the original 
ones. The role played by the conjugate field is clarified in 
this context, providing an alternative measurement foun- 
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dation of the formalism and results developed in Refs. 
[H,!!^. Moreover, some general results and conclusions 
follow from the measurement scheme. 

A general analysis is developed by assuming a renewal 
property flHij. Hence, the counting process is character- 
ized by a probability distribution, or waiting time distri- 
bution (WTD), that defines the statistics of the random 
time intervals between consecutive events. On the basis 
of the conditional measurement scheme, we show that the 
thermodynamic frame can also be related to an underly- 
ing renewal counting process whose WTD is the original 
one multiplied by an exponential decay scaled by the free 
energy function. Furthermore, it is demonstrated that 
scale invariant points, a shift closure property, intermit- 
tent and (thermodynamic) finite-size effects can be de- 
veloped under the renewal hypothesis. While no general 
results can be formulated for non-renewal counting pro- 
cess, we show that similar conclusions can be obtained for 
process defined by a stochastic WTD. This non- renewal 
case arises in different situations such as quantum optical 
and neuronal ones 
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systems [Ic 

The paper is outlined as follows. In Sec. II, after re- 
viewing the LD approach in the context of counting pro- 
cesses, we define and analyze the underlying statistical 
mechanics associated to the thermodynamic formalism. 
In Sec. Ill we study the central case of renewal counting 
process. In Sec. IV we analyze non-renewal processes 
defined by a stochastic WTD. In Sec. V we provide the 
conclusions. 



II. LARGE DEVIATION THEORY OF 
COUNTING PROCESSES 

A counting process can be statistically characterized by 
a set of probabilities {Pn{t)}^^Q, satisfying < P„(i) < 
1, and the normalization X]J^o-^"(^) ~ ^- Each P„(i) is 
the probability of occurrence of n-events up to time t. As 
in Refs. [1^, [2^ , here we are restricting our analysis to 
unidirectional counting processes, i.e., n > 0. 

While the short time behavior of the counting prob- 
abilities depends on each specific case, their asymptotic 
behavior may assume a universal structure. For example, 
CLT allows to approximate their long time regime with 
a Gaussian distribution 0, Q 



lim P„(i) 



clt 
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where n{t) is the average number 

oo 



(1) 



(2) 



n=0 



while Anf corresponds to the second cumulant 

oo 

A^, = J2{n - n{t)fPn{t) « {{An'))t. (3) 



((n)) and {{An?)) are the asymptotic growing rates. For 
a fixed time, the prediction of CLT is only valid up to a 
given accuracy in a central region of the counting prob- 
abilities. By assuming valid an exponential asymptotic 
behavior 



lim Pn{t) 

t—^OO 



exp[-i(^(-)], 



(4) 



LD method allows to describe the long time regime be- 
yond the Gaussian approximation. On the basis of a 
saddle-point approximation p^ . the function tp{n) can 
be obtained from a Legendre-Fenchel transformation 

(p{n) — max[8(s) — sn], (5) 

s 

together with the inversion formula 

0(s) = min[(^(n) -|- sn]. (6) 

n 

The function 0(s) is defined by the asymptotic behavior 
lim Z{t, s) w exp[-te(s)]. (7) 

i— f oo 

Here, s) is the characteristic function associated to 
the counting probabilities 



Z{t,s) = Y,Pr,{t)e 



(8) 



where s is a real dimensionless parameter. Hence, 6(0) — 
0. From relations ([5|) and ([6]) one deduces that any of the 
functions or 0(s) provide a complete description of 
the asymptotic counting statistics. 

CLT is covered by the LD approach. From Eq. ([J) the 
LD rate function ^p{n) can be approximated as 



ip{n) 



cu {n-{{n)))' 



2{{An^)) ■ 
Furthermore, Eq. ^ allows us to obtain 



clt 



e(s) « s((n))--((An^}). 



(9) 



(10) 



n=0 



This expression can be read as a Taylor expansion up to 
second order in s of 0(s), which in turn corresponds to a 
quadratic fitting of (p{n) around its minimal value, n = 
{{n)). Thus, the Gaussian prediction of CLT follows from 
a quadratic approximation of both ip{n) and 6(s). The 
LD method takes into account all higher contributions 
related to higher cumulants. 



Underlying thermodynamics and statistical 
mechanics 



Legendre-Fenchel transformations ([S|) and (jH]) suggest 
to analyze the LD method within an equilibrium ther- 
modynamic frame [1^ [2^. In fact, the functions (p{n) 
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and Q{s) can be read and related to the internal energy 
(U) and grand potential [6(s)] of a thermodynamic open 
system, which satisfy the relation [soj 

e{s) ~TS + sN^. (11) 

is the average "particle number," s plays the role of a 
(dimensionless) chemical potential (intensive parameter) , 
and the (dimensionless) temperature T can be taken as 
one. Here, we re-derive these relations and associations 
[2^ from an equivalent statistical mechanics approach. 
Hence, we search for a set of auxiliary (counting) proba- 
bilities {(i'n(i)}^o that maximize an entropy function S 
consistent with the thermodynamic relation ([Tl]) . 

As the LD approach applies in the asymptotic regime, 
we define a normalized entropy function as 

00 

5=- lim -^g„(01ogb„(t)]. (12) 

In order to find the probabilities {qn{t)}^=o, we maxi- 
mize the entropy function (at a fixed time) under the 
constraints of probability normalization X]^o (0 ~ 1 > 
fixed average (particle) number 

00 

N^^m)^ lini -Y,qn{t)n, (13) 

t-i-oo t — ' 

and (average) internal energy 

00 

V - ((£;)) ^ - lim - ^ g„(i) log[P„(t)]. (14) 

n=0 

By using the method of Lagrange multipliers [s^l, after 
defining the free energy or grand potential [see Eq. ([T])] 

e(s) = - lim ilog[Z(t,s)], (15) 

t— 5-00 t 

and using the thermodynamic relation Eq. (|11|) . we get 
the asymptotic behavior 

hm g„(i) « exp{-t[(^)s ^(^) - e(s)]}. (16) 

In this derivation, we assumed units of energy such that 
the Lagrange multiplier associated to the internal energy 
constraint (inverse of the temperature) takes the value 
one (T — > 1), and s corresponds to the Lagrange multi- 
plier associated to the constraint (|13|) . i.e., the dimension- 
less parameter s plays the role of a chemical potential. 

Eq. (fT6|) is satisfied if we define the probabilities 
qn{t, s) = qn{t) at any time as [see Eqs. (H)) and ^] 

g„(i,s)EE— l^P„(t)e-^", (17) 
Z{t,s) 

where the "partition function" Z(t, s) follows from Eq. 
(|S]). These probabilities define an extra counting process. 



which is parametrized by the dimensionless chemical po- 
tential s. Due to this dependence, the set of stochastic 
realizations consistent with {qn{t, s)}$^Lo is named as the 
s-ensemble [Ulii]. 

Relevant information about the fluctuations of the 
counting process defined by the set iPn(t)} ^ n is encoded 
in the thermodynamical response functions '30|, such as 
the first and second derivatives of the grand potential 
with respect to s 

m) = ^eis), ((A7V2)) = _|le(s). (18) 

OS os^ 

In fact, different (dynamical) phases, as well as thermo- 
dynamic transitions between them, can be established by 
analyzing the dependence of these objects with respect to 
the pseudo chemical potential [2^ [2^ . Due to this rea- 
son, the first response function is alternatively denoted 
as a dynamical order parameter" while s is quoted as 
its non-equilibrium conjugate (counting) field" [28| . 

Response functions (IT51) can be written as the normal- 
ized [limf_i.oo(l/t) • • •] average value and variance of the 
number of events associated to the set of probabilities 
{(7„(i, s)}^Q [see Eq. p^ ]. The same property is valid 
for higher objects. Hence, the s-ensemble provides an al- 
ternative basis for describing the thermodynamic frame. 
This fact defines the present approach. 

B. Realizations of the s-ensemble: conditional 
counting scheme 

The point process associated to probabilities 
{P„(i)}J^g is our observable system. Therefore, 
one has access to its ensemble of realizations. Never- 
theless, it is not clear which kind of point process leads 
to the counting process defined by the probabilities 
{qn{t, s)}J^g, Eq. pT)) . While its definition seems to be 
rather abstract, here we provide a measurement scheme 
for generating its ensemble of realizations from those 
of {Pn{t)}^^Q- The dependence of the thermodynamic 
frame on the conjugate field s is clarified by this result, 
which in turn answers the main issue raised up in the 
introduction. 

When s > 0, we realize that < exp[— s] < 1. There- 
fore, this factor can be read as a probability. In this 
case, the realizations associated to the counting proba- 
bilities {g„(i, s)}^o i'^ ti^6 interval (0, t) can be obtained 
as follows. Each realization corresponding to {Pn{t)}^^o 
is split in two channels (upper and lower arms, see 
Fig. 1) with probabilities exp[— s] and (1 — exp[— s]) re- 
spectively. If in the interval {0,t) all events are selected 
in the upper arm, that realization is taken as one of 
the set {qn{t,s)}^^Q [Fig. 1(a)]. If at least one event 
happens in the lower arm, the realization is discarded 
[Fig. 1(b)]. Therefore, each realization of {Pn{t)}^^o 
having n— events is selected as one of the {qn{t,s)}^^Q 
with probability exp[— sn] and discarded with probabil- 
ity ELi (fe)e-''("-'=Hl - e"")'' = (1 - e"*"). This con- 
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FIG. 1: Conditional counting scheme that leads to the realiza- 
tions of {gn(i, s)}^i!Lo- Each event of {-Pn(t)}^^Lo (left vertical 
lines) is selected with probabilities exp[— s] and (1 — exp[— s]) 
in the upper and lower arms respectively. In (a) all events 
went in the upper arm, leading to a s-realization. (b) Real- 
izations with at least one event in the lower arm are discarded. 

ditional selection corresponds the factor exp[— sn] in the 
definition ([T7|) . Furthermore, Z{t,s) can be read as the 
fraction of the total realizations of {Pn{t)}^^Q that are 
taken as realizations of the set {q-nit, s)}^^q. Therefore, 
it provides the right normalization associated to the con- 
ditional selection. Alternatively, it can be read as the 
probability of not occurring any event in the lower arm 
in the interval {0,t). 

When s < 0, it follows exp[— s] > 1. Therefore, the 
previous scheme does not apply. Nevertheless, a similar 
interpretation can be established. We realize that for any 
value of s it is possible to invert Eq. ((T7)) as 

Pn{t) = —^^qn{t,s)e^\ (19) 

where we have defined the partition function 

oo 

Z,(t,s)^^g„(t,s)e-^". (20) 

n=0 

Notice that the normalization of the counting probabili- 
ties implies the relation Zq{t, —s) = 1/Z{t, s). For s < 0, 
Eq. has the same structure than Eq. ([T7]). Hence, 
even in this case, the previous measurement scheme 
makes sense. Nevertheless, for s < the input signal 
must correspond to the realizations of {qn{t, s)}^o while 
the the output signal delivers the realizations defined by 
the set {Pn{t)}^=o- Although this association is exact 
(s < 0), it does not provide a procedure for generating 
the s— ensemble from the original counting process. In 
fact, it solves the inverse problem. 

The relation between the original ensemble and the s- 
ensemble can be extended to two arbitrary values, s and 
s', of the pseudo chemical potential. We write 

qn{t, s) = -J-—q^it, s')e-(^-^>, (21) 
Z{t, s, s') 

where the normalization here reads 

oo 

Z{t,s,s') = J2Qn{t,s')e-^'-''^^. (22) 

n=0 



For s > s' , the realizations of {qn{t,s)}^^Q follow 
from those of {<;n(i, sOl^o applying the conditional 
measurement procedure sketched in Fig. 1. The set 
{P„(i)}^Q is associated to s' = 0. Furthermore, it is 
simple to realize why it is not possible to build up a 
measurement scheme for obtaining the realizations of 
{(j„(t, s')};^o from those of {quit, s)}'^^q, i.e., the case 
s < s'. In fact, the realizations of the last counting pro- 
cess follow from a subset of the realizations of the former 
one, implying a lack of information which cannot be re- 
covered from any measurement scheme. 

From the definition of the conditional measurement 
scheme (Fig. 1) and the previous analysis, it is simple 
to realize that the thermodynamic response {{N)) [Eq. 
(|13p ] is a strictly decreasing function of s, implying 

lim {{N}) = 0, lim {{N)) = oo. (23) 

s— ^-t-oo s— > — oo 

These constraints imply the concavity of the free energy 
function 6(s) and the limits 

lim e(s) = |C|, lim e(.s) = -oo, (24) 

s— ^-|-oo s— > — oo 

where C is a real constant. Consistently with the behav- 
ior of {{N)), in general {{AN'^)) is neither an increasing 
nor a decreasing function of s. 

The function Z{t, s) measures the fraction of the num- 
ber of realizations from the original counting process 
that are taken for the s-ensemble. The asymptotic be- 
havior ([7]) implies that a reliable measurement of the 
set {gn(i, s)}^o [^1- C3] ^^'^ response functions 
involves a number of realizations of {P„(t, s)}5^q that 
scales as exp[t*0(s)], where t* is larger than the transi- 
tory time interval. From a numerical point of view, this 
condition is very restrictive. In fact, it can take a long 
time before a realization is generated. In spite of this 
severe limitation, below (Fig. 3) we show that the mea- 
surement scheme can in fact be numerically implemented 
for moderate values of s, i.e., s « 0, such that the product 
t*Q{s) does not assume very high values. 



III. RENEWAL COUNTING PROCESSES 

While the thermodynamic potential Q{s) provides a 
complete description of LD fluctuations, its functional 
form depends on each specific case. In contrast, some 
general results can be established from the statistical me- 
chanics formulation based on the set {qn{t, s)}^^q. Here, 
we use this advantage for providing a general character- 
ization of an arbitrary renewal counting process 0, Q . 

A renewal point process is defined by a WTD satisfying 
w{t) > and the normalization dt'w{t') = 1. It de- 
fines the probability distribution of the (random) time in- 
terval between consecutive events. Therefore, the count- 
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ing probabilities {P„(t, s)}5^q can be written as 31| 



PoW = 1- / dt'w{t'), (25a) 
Jo 

Pn{t) = / dt'w{t~t')Pn-l{t'). (25b) 

Jo 

By working in a Laplace domain, it is simple to derive 
the equivalent evolution equations 



-Po(i) = - dt'K{t-t')Po 



(26a) 



dt 



Pu{t) 



dt'K{t - t')[P^{t') - Pn-l{t')]. (26b) 







The memory kernel K{t) is defined in the Laplace do- 
main, g{u) = dt(7(i) exp[— ut], as 



K{u) = 



w{u) uw{u) 



Po{u) l-wiu)' 



(27) 



While the kernel K{t) does not have a straightforward 
physical interpretation, its time integral does. After 
defining f{t) = dt'K{t'), or equivalently in the Laplace 
domain f{u) = K(u)/u, from Eq. (P7|) we get 



fiu) 



w(u) 
1 — w(u) ^ 



(28) 



By writing this expression in the time domain in terms 
of successive convolutions of w{t), we realize that f{t)dt 
provides the probability of observing an event in (t, t+dt) 
independently of the occurrence of any extra events in 
previous time interval {0,t) [3l|. In contrast, w{t)dt has 
the same interpretation under the condition that not any 
event occurs in (0,t). 



A. s-ensemble 

As the renewal property implies the absence of mem- 
ory between consecutive events, the s-ensemble p7)) can 
easily be characterized. From Eq. the evolution of 

the partition function ([5]) reads 



— Z(i,.s) = -(l-e-") / dt'K(t~t')Z(t',s), (29) 
dt Jq 



whose solution in the Laplace domain is 
Z{u, s) 



u + K{u){l - e-") 



(30) 



By writing the survival probability in the Laplace do- 
main as Pq (u) = [u + K{u)]-'^ [see Eq. it follows 
that Z(u, s) with s > can be read as the Laplace trans- 
form of a survival probability defined with the normal- 
ized kernel K{u) — (1 — e~^)K{u), which in turn implies 



f{u) — ^ (1 — e~^)f{u). Therefore, we confirm that in fact 
Z{t, s) can be identified with the survival probability of 
the lower arm of the measurement scheme of Fig. 1. 

From its definition and using the renewal property, 
Eq. ([25)1 . Z{t, s) can alternatively be written as 

Z{t,s)=PQ{t)+ [ dt'w{t-t')e-^Z{t',s). (31) 
Jq 

Then, as qQ{t,s) — Pa{t)/Z{t,s) it follows 



9o(i,s) = 1 - / dt'wq{t,t',s). 



(32) 



Furthermore, from Eqs. (IT71) and (P5|) . for n > 1 we get 



qn{t,s)= / dt'wq{t,t',s)qn-i{t' ,s), 
Jo 



where we have defined 



(33) 



(34) 



This expression, together with relations (|32]) and (|33l) . 
imply that the s-ensemble is a non-stationary renewal 
counting process [compare with Eqs. (f25|) ]. In fact, 
here the WTD Wq{t,t',s) not only depends on the in- 
terval between successive events, but also on the time 
of the last event. Consistently with the measurement 
scheme of Fig. 1, this dependence is introduced by 
the partition function factors. By introducing the in- 
terval T = t — t' between two consecutive events, we 
write 'Wq{t,t — r, s) = e~''w{T)Z(t — t, s)/Z(t, s). In 
this way, Wq{t,t — r, s) can be read as the probabil- 
ity density for observing an event at t given that the 
last one occurred at time (i — r). The normalization 
limt_j.oo dTWq{t,t — T,s) — 1, follows from Eqs. (PT|) 
and ([5^ after using that limt^^oo qo{t, s) = 0. In the long 
time regime, where the partition function can be approx- 
imated by an exponential decay, Eq. ([7]), we get 



W?°{t,s) = Ihxi Wq{t,t — T, s), 

= e~*i(;(T) exp[T0(s)] 



(35) 



Therefore, in the regime where the LD approach applies, 
the counting process defined by the set {qn{t, s)}5^o be- 
comes a stationary renewal one, which in turn is con- 
trolled by the WTD 'w^{t, s). This is one of the central 
results of this section. Notice that it is valid for any 
value of s. On the other hand, one can also associate to 
the s-ensemble a non-conditional probability distribution 

fq{t, t', s)^Wq{t, t', S)+J^ dtlWq{t, tl , s) fq{tl , t' , s), (36) 

which defines the probability density of observing an 
event in (t, t + dt) given that one occurred at t' and in- 
dependently of the occurrence of any extra events in the 



6 



time interval {t,t'). By solving Eq. (p6)) iteratively, after 
some calculations we get 

f^{t,t',s)^e-^r{t~t',s)^l^, (37) 

where F{t, s) is defined in the Laplace domain as 

K(u) 



F {u, s) — K(u)Z{u, s) 



u + K{u){l ~ e'^) 



(38) 



Consistently, fq{t,t',0) = f{t — t'). Furthermore, in 
correspondence with Eq. (1551) . we define f^{T,s) — 
limj_j.oo fq{t, t ^ T,s), which reads 



(39) 



From its definition 31], it follows the relation 
\imr^^f^{T,s) = ((TV)), Eq. (USD. 



B. Free energy and response functions 

The underlying renewal property of the s-ensemble al- 
lows us to get some general expressions for the free en- 
ergy and its response functions. The Laplace transform 
of w^{t,s) [Eq. ([55)) ] with respect to r can be written 
as w^(w, s) = e~^w{u — 0(s)). As this function satisfies 
the normalization dTw"^ {t, s) = 1, or equivalently 
w'^{u,s)\u=o = 1, we obtain 



w{u - e(s))|„=o = 



(40) 



By using this relation and the asymptotic values (1^^ one 
can deduce the relation 



lim e(s) «7o(l-e-"«^). 



(41) 



valid when w{t) has a characteristic short time scale, i.e., 
limu^^ w{u) ~ (7o/w)i/"°. 

The Laplace transform of w{t) also provides an equa- 
tion for getting the free energy function. In fact, it is the 
smaller solution of Eq. This relation is valid for any 

(ergodic) renewal process. Notice that the same kind of 
relations follow from definition ([T])- In fact, by using the 
residues theorem, 0{s) can be determined as the smaller 
root of the denominator of Eq. ([50)1 . 

As the process defined by the set {qn{t, s)}'^^q is a 
renewal one in the stationary regime, the thermody- 
namic response functions [Eq. dlSp ] can be derived from 
?«^(t, s). In fact, by defining the time interval moments 



(42) 



T„(S) = / W^iT,s)T"dT, 

Jo 

and by using the renewal property, it follows 
1 



m) 



im^))-'-^^^^, (43) 



where we have used that {{N)) and {{AN'^)) measure the 
linear growth of the first two cumulants of the counting 
process defined by the set {g„(t, s)}5^o (see page 496 of 
Ref. "sill). Similar relations are valid for higher response 
functions, i.e., u;^(r, s) completely determine the ther- 
modynamic formalism. 



C. Applications 

Below, we show some applications that rely on the pre- 
viously developed results. 



1. Scale invariant renewal process 

In Ref. [2^, , a scale invariant property was found 
in a class of photon counting process. It is defined by the 
condition 



m) 



(44) 



where the response functions follow from Eq. (jl8l) . and 
a is a positive real constant, < a < 00. Hence, the 
normalized fluctuations do not depend on s. The count- 
ing processes {P„(i)}Jf^o ^nd {g„(t, s)}Jf=0' ^^^y value 
of s, are then sub-Poissonian for < a < 1, while for 
1 < a < cx), becomes super-Poissonian 0-11] ■ Here, we 
demonstrate that the scale invariant property is always 
satisfied by a class of (gamma) renewal process [3-3 j 
which are defined by the WTD 



w{u) - 

In the time domain it reads 
wit) - — ^ 

r(i/a) 



l/a 



(45) 



(46) 



The average interval between consecutive events is given 
by 'w{t)tdt = (aj)~^ . By working the expressions 
P5|) in the Laplace domain, the counting probabilities 
read 



Pn{t) = 



r(^,70 r[g,7^] 
r(^) r(^) 



(1 - 6no), (47) 



where r[fc, x] is the incomplete gamma function r[fc, x] = 
z^'^^e^^dz, and T{k) is the Euler gamma function, 
i.e., T{k) — T[k, 0]. In the long time regime, ^ 1, it is 
valid the approximation 



Pn{t) 



1 



-ilt)'- 



-exp[-7t]. 



(48) 



^ a 

From this expression, it follows the LD function [Eq. (U) 

n r / n 

— 1 - log — 

q;7 L \Q!7^ 



ip[n) = 7|l 



(49) 
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{n/t) n. Its Legendre-Fenchel transformation [Eq. ([6])] 
leads to the grand potential 



(50) 



Consistently, this result straightforwardly follows from 
Eq. (gOl)- Alternatively, Eq. ^U\j leads to the partition 
function 



1 (^ + 7)V"-y/" 
Ziu, s) = —- —. —. 



(51) 



From the residues theorem, the free energy ()50p is re- 
covered. In fact, u = —Q{s) cancels the denominator 
of this expression. The grand potential ([50]) implies the 
thermodynamic response functions [Eq. ([T8|] 



((AiV^^) = a2^e-"^ (52) 



Trivially, these expressions satisfy the invariant property 
Eq. (pj). 

While it is not possible to find a simple expression 
for the non-stationary WTD Wq(t,t',s) [Eq. ([M])], the 
stationary one [Eq. ([35)) ] reads 



<(r,s) 



1 



ni/a) 

where the renormalized rate reads 



7sexp[-7sTj 



7s 



7exp[— asj 



(53) 



(54) 



The stationary WTD w^{t,s) is the same as the orig- 
inal one [see Eq. (j46p ] under the replacement 7 ^■ 
7e~"*. Thus, the counting process defined by the set 
{quit, s)}5^oi asymptotic regime, can be obtained 

from the original one after a time rescaling 



(55) 



lim qn{t,s) ~ lim P„{te'°"). 



This condition was derived from the WTD Eq. (|45|) . 
Nevertheless, it is simple to prove that its fulfilment is 
sufficient to guaranty the invariance property (|44|) . By 
using that Zg{t, -s) = 1/Z{t, s) [see Eq. (gHl)], from ^ 
it follows the equivalent condition 



e(-s) = -e(s)e° 



(56) 



which in turn allows to write an equation for {d/ds)0{s). 
Under the constraints ([55)) . continuous derivatives in s = 
0, and S(0) = 0, it has a unique non-null solution that is 
given by Eq. ([50]) . Notice that this derivation does not 
rely on the renewal property. Nevertheless, the previous 
one guarantees that {-Pn (i)}n^o asymptotically converges 
to a renewal process defined by the WTD ps)) . 

When a = 1, the WTD becomes an exponential 
one, w{t) = 7exp(-7t). Thus, f{t) = -f9{t) [Eq. p8|] 
where 9{t) is the step function, implying a local in time 
evolution [K{t) = ^5{t)] of the counting probabilities [Eq. 
CT ]. Their solution read P„(t) = (7^)" exp[-7i]/n!, i.e.. 



the very well-known Poisson process [ll-Ql . The partition 
function [Eq. (|5])] reads 



Z(t,s) = exp[-7t(l - e^")] 



1, 



(57) 



which in turn implies the free energy ()50p and its asso- 
ciated response functions with a = 1. The s— ensemble 
[Eq. (IT7|) ] at "any particular time" can easily be writ- 
ten as qn{t,s) = {e~'^jt)^exp[—e~^jt]/n\. Consistently, 
these probabilities also correspond to a Poisson process 
with the normalized rate 7 — > e^'^j. Its associated non- 
stationary WTD [Eq. (p4)) ] is stationary at all times 

Wg{t,t',s)^{-fe-')cxp[-{-fe-'){t-t')], a = 1, (58) 

i.e., at any particular time it only depends on the differ- 
ence {t — t'). This property is only valid for the Poisson 
process, a — 1. 

Scale invariant photon-emission process. In Ref. (28j 
it was analyzed the case of a two-level fiuorescent system, 
where an external laser field induces a continuous emis- 
sion of photons at random times. The photon counting 
process is a renewal one, being defined by the WTD 



w{u) 



7/2 



-7/2' 



- W7 - 



(59) 



Here, u is the Laplace variable, 7 denotes the natural 
decay rate of the system while the Rabi frequency fl 
measures the system-laser coupling. The average time 
is J^w{t)tdt = [-/n^/if + 2n^)]-\ When 217 = 7, 
it was found the invariant scale property (|44p with 
{{AN^}) / {{N}} = 1/3. In this situation, §^ becomes 



w{u) = 



n 



(60) 



which in the time domain leads to w{t) = 
ijl'^t^ exp[— fit]. These expressions and the previ- 
ous results confirm that in fact a = 1/3 [see Eq. (|45|) ]. 
showing the consistence of the present approach. 



2. Shift closure property 

With {f} we denote the characteristic parameters 
(rates, etc.) that determine a given counting process 
{P„(t)}^g. Evidently, its associated free energy depends 
on {^}, property denoted as [6(s)]{^}. We define a shift 
closure property by the condition 

[0(s)]{?(.o)} = [e(s + so)-e(so)]{c}, (61) 

where {^(so)} are new characteristic parameters that de- 
pend on So. Note that [9(s + so) — 6(so)]{^} is the free 
energy function of the counting process {gn(i, so)}5^o- 
In fact, all response functions, i.e., its derivatives with 
respect to s [Eq. p^ ]. are the same as those of 0(s) 
shifted by — so. Hence, the condition (pT|) tells us that 
{qn{t, 50)1^0 (tli^ so-ensemble) belongs to the family of 
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processes obtained from {Pn{'t)}'^=o by a change of char- 
acteristic parameters, {£,} {C(so)}- Consequently, the 
fulfilment of Eq. (|6T|) allows us to generate the realiza- 
tions of the s-enscmble by changing the parameters of 
the original counting process. Equivalently, it implies 
that the change of parameters {^} — > {'f(so)} in the def- 
inition of {Pn{t)}^^Q only produce a shifting of all of its 
response functions. 

Scale invariant processes, Eq. ([50)) . satisfy the shift clo- 
sure condition with the mapping {7} — >■ {76""*"} [Eq. 
([M])]. In general, it is very difficult to check the validity 
of (pT|) . Nevertheless, for renewal processes it can alter- 
natively be written in terms of the stationary WTD ((35)) 
as 

<(r,s) = [«;(r)]{4(,)}. (62) 

Hence, w'^{t, s) follows from w{t) after replacing {^} 

(s)}. When satisfied, this equation is not only equiva- 
lent to (j6ip but also provides a simple way for obtaining 
the parameters mapping {f} — {f(s = so)}- 

A broad class of WTDs satisfy the shift closure prop- 
erty. For example, consider in the Laplace domain 
w{u) — Y[i Ai/(w + Xi). By using the normalization con- 
dition (|40|) . we get 

w^{u,s)=l[ y-') A,(s) = A,-e(.s). (63) 

Notice that w^{u, s) have the same structure than w{u) 
with the rescaled parameters {^(s)} = {Ai(s)}. This sim- 
ple result guarantees the validity of Eq. (pTj) . In gen- 
eral, Aj may be complex numbers (i?e[Ai] > 0), which 
in turn are functions of experimental parameters. Eq. 
(|59l) has this structure, with Xi — Xi{j, fJ^). In this case, 
we checked that w'^{u,s) can be written as wlu) af- 
ter replacing 7 ^ 7(5) and fl^ i^^(s), where 7(5) — 
7-2e(s), and n'^{s) = n^ + e(s)[e{s)--f], where Q{s) is 
defined by Eq. (35) of Ref. Notice that the param- 
eter transformation satisfies 7^(s) — 4r2^(s) = 7^ — 4i7^. 

Another class of WTDs that satisfy the shifting condi- 
tion ([62]) is given by w{t) = J^kPklke''^''*, with J2kPk = 
1, < Pk < 1. Thus, in this case {£_} = {jk,Pk}- From 
Eq. dSSI), we get 

%°°(^.^) = (64) 

k 

where the parameters {C('S)} — {'Jk{s),Pk{s)} are 

7fc(s) = 7fe - ©(•«), Pkis) = Pk ^ '1) (65) 

7fc - e(s) 

with 7fc(s) > and J2kPk{^) = 1- Notice that w'^{t,s) 
can also be written as a superposition of Poisson distri- 
butions, which guarantees the fulfilment of (|61l) with the 
mapping {-fk,Pk} {7fc(so),Pfc(so)}, Eq. ([65]). 

In order to exemplify the previous results, we consider 
a bi-exponential WTD 

w{t) = (1 - p)7>e~''>* + p7<e-''<*, (66) 




FIG. 2: Free energy function (|69|) (full line) in terms of s. 
The parameters of the WTD [Eq. ^] are 7</7> = 0.007 
and p = 0.015 [7/7> = 0.9815]. The dotted line corresponds 
to Eq. (|61|) with so = —0.1. The short-dashed line is the 
estimation based on the CLT, Eq. (|10|l . 

where 7> and 7< are the characteristic rates and the 
parameter p satisfies < p < 1. The kernel (|27l) defines 
the evolution of the counting probabilities {Pn{t)}^^Q 
[Eq. (HH)]. It reads 

K{t)=j{Sit)-PeM-Vt]}, (67) 

where the characteristic parameters are 

7 = (1 -_p)7> +P7<, (68a) 
/3 = {(l-p)7>+P7<-7'}/7, (68b) 
^7 = (1 -Ph< +P1>- (68c) 

The grand potential [Eq. (fT5|) ] follows from the relation 
(|40p. Equivalently, by using the residues theorem it can 
be obtained from the smaller root of the denominator of 
Eq. dSOI), + m[7(1 - e"") +ri\ + {f]~ /3)7(1 - e"^) = 0. 
We get 

e(s) = 2(l-e-) + |-i{[,7-7(l-e-^)p 

+4^7(1 -e"")}^^^ (69) 

In Fig. 2 we plot S(s) (full line) as a function of 
the chemical potential s for a particular set of param- 
eter values of the WTD jM])- The limits ^ are satis- 
fied. The dotted line corresponds to [9(s -f so) — ©(so)], 
with So = —0.1. Consistently, we checked that the same 
(shifted) free energy function follows from Eq. ((69)) 
under parameter replacements {^} = {p, 7>,7<} 
{p(so),7>(so),7<(so)} = {C(so)}, defined by Eq. ([65]). 
This example explicitly shows the meaning of the shift 
closure condition. 

In Fig. 3 we plot the response functions, Eq. ((T8|) . 
Consistently, they fulfill Eq. ((15|) . We also checked 
that the shifted response functions coincide with the 



9 



1.5- 



1.0' 



0.5- 



0.0 




(a) 



exp[-s] 



-0.2 0.0 0.2 
S 

— 1 1 , 1 ' I 

■0.4 -0.2 0.0 0.2 0.4 




-1 — I — 1 — I — 1 — I — 1 — I — 1 — I — I- 

-0.1 0.0 0.1 0.2 0.3 0.4 0.5 




2000 3000 



FIG. 3: (a) First and second (inset) response functions, Eq. 
(|18p . The parameters of the WTD are the same than in Fig. 
2. The dotted hnes correspond to a set of different parameters 
values obtained from Eq. (|65[l with sq = —0.1 and so = —0.2. 
The circles correspond to the numerical values obtained from 
the measurement scheme of Fig. 1 with the transformations 
s = -> s = 0.001, 0.015, 0.05. The short-dashed hnes are 
the estimations based on the CLT, Eq. (|10|) . The long-dashed 
lines correspond to the response functions of a Poisson pro- 
cess with rate 7 (see text), (b) Normalized second response 
function {{AN^)) / {{N)) . 



derivative with respect to s of [8(s)]{^(5q)}. In Fig. 3(a) 
all plotted functions are normalized by the rate 7, Eq. 
(|68ll . which is different for each set of parameter values, 
s — 0, —0.1, —0.2. Hence, the shifted response functions 
differ from the original one (s = 0) in a multiplicative 
factor. This effect is absent in the normalized second 
response function {{AN^}} / {{N)) , Fig. 3(b). 



3. Intermittent renewal process and finite-size effects 

In order to enlighten the LD method, in Figs. 2 
and 3 we show the predictions corresponding to the 

CLT (short-dashed lines), i.e., e(s) = s{{{N))\,=o) - 
{s^/2){{{AN'^))\s=o), Eq. ([10]). Thus, the first response 
function has a linear dependence on s and the second 
one becomes a constant. While the first approximation 
is valid in a small region around s w 0, the second one 



FIG. 4: Stationary survival probability (|70p as a function of 
time for different values of s. The inset corresponds to a log- 
log plot. The intensity realizations (see text) ploted in (a), 
(b), and (c) correspond to s = +0.15, s = 0., and s — —0.4 
respectively. The parameters of the WTD (|66[l are the same 
than in Fig. 2. 



only applies at s — 0. 

The previous analysis implies that the CLT approxi- 
mation could miss a phase transition within the thermo- 
dynamic approach. In contrast, the response functions 
plotted in Fig. 3 develops properties consistent with 
finite-size effects in a first order transition [H, H^l- In 
fact, the first response function changes abruptly around 
s « 0, while the second one develops a narrow peak. As 
shown in Refs. [2^, [2^ these features can be related to an 
intermittence (blinking) property of the counting process. 
Here, we show that renewal processes, independently of 
the underlying dynamics, may also develop similar phe- 
nomena. 

The main features of the response functions shown in 
Fig. 3 can be easily interpreted in terms of the s-ensemble 
dynamics, which in turn is defined by the stationary 
WTD ([35t . In Fig. 4, in order to show the s-dependence 
of w'^{t, s), we plot its associated survival probability 



Jo 



(70) 



or equivalently w^(t, s) = —{d/dT)q^{T,s). This func- 
tion allows to generate the ensemble of realizations con- 
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sistent with the asymptotic behavior of the counting 
probabihties {qn{t, 3)}"^^^. In fact, the time interval be- 
tween consecutive events can be determined by solving 
the equation q^(T,s) — r, where r is a random number 
in the interval (0, 1). 

The previous algorithm allows us to check the mea- 
surement scheme defined in Fig. 1. We numerically im- 
plemented the conditional selection to the generated re- 
alizations (s — 0) and calculated numerically {{N)) from 
the asymptotic behavior ([T3|. ((AiV^)) can be obtained 
in a similar way as the growing rate of the quadratic fluc- 
tuations. In Fig. 3 we show the results (circles) for the 
transformations s = — >■ s = 0.001, 0.015, 0.05. Even for 
values of s where the predictions of the CLT do not ap- 
ply, the numerical estimations agree with the theoretical 
results. 

In Fig. 4, for each chosen value of s, we also plot- 
ted the intensity I{t) associated to an arbitrary stochas- 
tic realization of the s-ensemble. They were obtained 
numerically from q^{T,s). The intensity is defined as 



(t)]/A, where 

{t) is the stochastic number of events up to time t and 
A define an adequate discretization time. Consistently, 
one can always find a A where I{t) fluctuates around 
{{N)) (Fig. 4(a) and (c), horizontal dotted lines). Nev- 
ertheless, for s — 0, i.e., for the counting probabilities 
{P„(i)}5^Q, one can also find a discretization time where 
the intensity develops an intermittence phenomenon, i.e., 
there exist successive time intervals where either many 
events happen or none happen at all. 

In the thermodynamic picture, the active and inac- 
tive periods are associated to two different phases, where 
finite-size effects lead to the transitions between them 
[29| . In the present approach, the intermittence phe- 
nomenon can be related to the WTD w{t) [Eq. 
For the chosen parameter values, which satisfy p ^ 1 
and 7</7> ^ 1, 'w{t) develops two different time scales, 
where the weight of one of them is much higher than the 
other one. This feature is evident in Fig. 4 [s = 0). 
Hence, the origin of the intermittence property follows 
straightforwardly from the previous defined algorithm 
11| . The active and inactive time periods arise respec- 
tively from the short and long time behaviors of the 
WTD. Consistently, the intensity of the active periods 
fluctuates around 7> « 7 (dotted line in Fig. 4b). 

The intensity realizations show that when s departs 
from zero, the intermittence phenomenon is lost. In the 
thermodynamic picture this property means that one of 
the two phases becomes dominant or stable. Specifically, 
the active and inactive phases are dominant for s < and 
s > respectively [see Fig. 3(a)]. Here, these conclusions 
can be related to the structure of the stationary WTD 



w^{t,s) [Eq. dai])]. 

The free energy G(s) goes to higher negative values for 
increasing negative values of s [see Fig. 2]. Hence, for 
s < the two time scales property is lost, and the ex- 
ponential term exp[— r|0(s)|] becomes dominant in Eq. 
In consequence the s-ensemble, by using Eq. (1411) . 



can be approximated by a Poisson process ^^(t, s) « 
7exp[— 7r], i.e., the active phase. The Poisson property 
explains why for s < all (normalized) response func- 
tions, independently of sq, converge to the same asymp- 
totic curve [Fig. 3(a)], {{N)) /7 « ((A7V2^)/7 « exp[-s]. 

For increasing values of s > the free energy is posi- 
tive, inducing a slower decay of w'^{t, s). As can be seen 
in Fig. 4, the survival probability maintains two different 
time scales, nevertheless their relative weights are similar. 
Therefore, in this case the intermittence phenomenon is 
also lost. Furthermore, as the weight of the larger time 
scale increases with s, the intensity is diminished, i.e., 
the inactive phase becomes dominant for increasing s. 

From the previous analysis, we conclude that the 
abrupt changes of {{N)) and {{AN^)) around s « 
refiect both the loss of the intermittence property for 
\s\ > 0, and the transition between the active and in- 
active phases. In the thermodynamic approach, the 
behavior around s « is read and estimated from a 
thermodynamic- finite-size analysis [1^. Those results 
can be straightforwardly applied in the present context 
after decomposing the partition function, Eq. 1^^, as 
Z{t, s) = ZA{t) + Zi{t), where the evolution of each con- 
tribution (Active, Inactive) can be written as 



dZA{t, s) 
Jt 

dZi{t,s) 
Jt 



-7(l-e-^)Z^(i,s) 

-TAZA{t.s) + TiZi{t,s), (71) 

+TAZA{t,s)-TiZi{t,s). 



Here, F^ = j3, Ti ~ rj — (3 = 7>7</7, and the initial 
conditions must to be Za{0,s) — 1, Z/(0,s) — 0. This 
splitting allows to read the counting process {-Pn(i)}^o 
as a stochastic modulated Poissonian one [see Eq. (|57|) ] , 
whose rate 7st(i) at random times adopts the values 
jA = J and 7/ = 0. The probabilities of assuming each 
value are governed by a classical master equation with 
transition rates F^ and P/. 

In a slow modulation limit [1^, i.e., when Ta/'^ ^ 1 
and P//7 ^ 1, Eq. (|71l) leads to an intermittence phe- 
nomenon. In fact, these conditions guarantee the exis- 
tence of periods of time with many consecutive events 
and periods where no events happen at all. The pre- 
vious inequalities can be satisfied by demanding p <C 1 
and 7< ^ 7>, such that 7</7> ~ 0{p). Hence, from 
Eq. (1551) it follows Pa ~ P7<, and P/ « 7<. The pa- 
rameters chosen in Fig. 2 satisfy these constraints. On 
the other hand, from Ref. [1^ we know that the transi- 
tion rates control the thermodynamic finite-size effects. 
In fact, the wide Up of the peak of the second response 
function {{AN'^)) [inset of Fig. 3(a)] here can be esti- 
mated as CTp ~ 2{rA -|- P/)/7 ~ 2p (see Eq. (66) in Ref. 
[2^). Therefore, in the limit p — 0, a first order tran- 
sition is approached, i.e., the first and second response 
functions converge to a discontinuous and a delta Dirac 
function respectively. Due to the validity of the shift 
closure condition, after a change of characteristic param- 
eters, Eq. (j65p . the transition may be observed at an 
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arbitrary value of the parameter s (see Fig. 3). 

In general, when the WTD involves more terms [Eq. 
([M)) ]. it is not possible to find an equivalent description 
such as the one given by Eq. (|71l) . Nevertheless, a similar 
approximated description can be formulated when the 
intermittence phenomenon develops. 

IV. NON-RENEWAL COUNTING PROCESSES 

While the renewal processes admit a relatively simple 
description, a general formalism for dealing with non- 
renewal counting process does not exist. In fact, the 
presence of memory between successive events does not 
admit a general description. Nevertheless, many cases 
can be mapped with a situation where the probability 
distribution for the interval between consecutive events 
(WTD) changes randomly after each event [l^ [l3l ■ 

We introduce a WTD w{T,Xb\xa) that, in contrast to 
the renewal case, also depends on an extra unobserv- 
able (hidden) variable x whose states are {xi}. With 
Xa and Xb we denote the states of the hidden variable 
before and after an event. Then, w{T,Xb\xa) not only 
provides the statistics of the time intervals but also de- 
fines the probability transition between the states of 
X. The most simple case corresponds to u;(t, x^jxa) — 
w{t, Xa)T{xb\xa), whcrc w(t, Xa) is a WTD parametrized 
by Xa and T{xb\xa) is the transition matrix between 
the hidden states. Processes with memory of previous 
events, processes with adaptation and double stochastic 
processes can be covered with this approach [l7|. Fur- 
thermore, a stochastic WTD may arise when describing 
quantum dissipative dynamics driven by classical fluctu- 
ations. In this last and Xb correspond to prob- 
abilities Pa and Pb of the initial and posterior hidden 
configurations [l3{ . 

By introducing a vector space spanned by the states of 
X, {|x)}, and by using the definition of w{T,Xb\xa), the 
counting probabilities can be written as 

Pn{t) = {l\P„{t)\p^n). (72) 

Vector \pin) denotes the initial probabilities of the hid- 
den variables, while (1| = (1, • • • , 1). The operators P„(t) 
follow from a convolution 

Pn{t)^ f din f"dtn-l--- f ^ dtiP„[t, {U}^^^], (73) 

Jo Jo Jo 
where the "joint probability operator" reads 

Pn[t, {t»}r=l] = Poit - tn)W{tr, - t„_i) • • • W (h) . (74) 

The survival operator is Po{t) = 1 — dt'W{t'), and 
W{t) has components 

{xj\W{t)\xk) = w(t,Xj\xk). (75) 

With the previous prescriptions, it is simple to derive 
similar expressions for the s-ensemble [Eq. (|17p ] 

q„{t,s) ^ {l\qn{t,s)\p^n), (76) 



where now the operators read 

Qn{t,s) = / dtn / dtn-1--- / [i, {^j} "Li , s] , 

Jo JO JO 

with the joint-probability operator 

Qn[t, {U}"^i,s] = qoit, tn,s)Wq{t, t^-l, s) ■ ■ ■ Wq{tl, , s) . 

The survival operator reads qo{t,t',s) — Po(i ~ 
t')Z{t' , s)/Z{t, s), while the waiting time operator is 

Wq{t,t',s)^e-^Wit-t')^^, (77) 

The partition function reads 

Z{t,s)^{l\Z{t,s)\p,n), (78) 
where its associated operator is 

oo 

Zit,s) = J2Pnit)e~'"- (79) 

n=0 

These expressions imply that the counting process de- 
fined by the probabilities {9n(i, s)}5^i has a similar 
structure to that of the set {P„(i)}^]^. In fact, the trans- 
formation of the waiting time operator [Eq. (j77p ] is sim- 
ilar to that of the renewal case [Eq. ([M]) ]. Hence, in the 
long time regime, it is also possible to define a station- 
ary waiting time operator W^(r, s) — limt^oo Wl^{t,t~' 
T, s), which reads 

W;^{t,s) = e-"W{T)exp[Te{s)]. (80) 

Thus, the non-renewal structure, in spite of the exponen- 
tial decay measured by 6(s), remains the same. Similar 
results to that obtained for renewal processes can be es- 
tablished from this last equation. 

V. SUMMARY AND CONCLUSIONS 

The LD method allows to characterize the statistics 
of a point process with a thermodynamic frame that is 
defined by the scaling rates of the counting probabilities 
and its characteristic function. The basis of the present 
analysis consists in finding an auxiliary process that max- 
imizes the entropy function of the thermodynamic ap- 
proach, Eq. (|17p . The conditional measurement scheme 
defined in Fig. 1 provides its ensemble of stochastic real- 
izations from those of the original one. This result gives 
an alternative measurement interpretation of the ther- 
modynamic frame and its conjugate counting field. Fur- 
thermore, it permits to relate the large fiuctuations of a 
counting process with the statistical properties of a sub- 
set of its realizations. 

Diverse general results were established under a re- 
newal property, i.e., when the successive events do not 
develop any memory between them. We established that 
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the s-ensemble, independently of the structure of the 
WTD of the counting process, is also a renewal process. 
An exponential decay measured by the free energy func- 
tion determines its (stationary) WTD [Eq. ([55]) ]. This 
function completely determines the thermodynamic for- 
malism. In fact, a general equation for the free energy 
can be established in terms of its Laplace transform, Eq. 
(HO)) . Similarly, the response functions follow from its 
normalized moments, Eq. (|43p . 

Different phenomena that appear in the thermody- 
namic formalism were analyzed and derived in terms of 
the developed approach. We demonstrated that an in- 
variance scale property, Eq. (|44)) . is satisfied whenever 
the s-ensemble, in the long time regime, follows from a 
time rescaling of the original one, Eq. (|55p . Equivalently, 
the invariance scale property arises when the counting 
process can asymptotically be approximated by a renewal 
one defined by a Gamma distribution, Eq. (HSI) . A shift 
closure property was introduced [Eq. (|5T|) ]. Its fulfill- 
ment implies that all possible statistical behaviors of the 
s-ensemble can be achieved and generated from the orig- 
inal one after an adequate redefinition of its characteris- 
tic parameters. Alternatively, it means that there exists 
a transformation of the characteristic parameters of the 
counting process whose sole effect is to shift all thermo- 
dynamic response functions in the s-direction. Different 
families of WTD guarantee the fulfillment of this prop- 
erty, Eqs. ((63)) and ((64)) . Finally, we showed that renewal 



processes can in fact develop intermittent phenomena, 
where the realizations of the process and the s-ensemble 
may be characterized by intervals of high and null count- 
ing rate (Fig. 4). The properties of the response func- 
tions, which approach a first order transition, were re- 
lated to the stochastic dynamics of the s-ensemble. While 
no general results can be established for non-renewal pro- 
cesses, we showed that similar considerations can be de- 
rived when the process is defined by a stochastic WTD, 
Eq. ((80t. 

The present approach gives a solid basis for under- 
standing the LD thermodynamic frame of a counting 
process. Its main advantage is the possibility of getting 
some general conclusions that apply to a wide class of 
point processes, independently of its specific structure. 
Extension to bidirectional counting processes, as well as 
non-ergodic ones may, in principle, be done along similar 
lines. 
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